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1. INTRODUCTION 

During processing of signal, unwanted parts of the signal are removed by use of a filter. This paper 
looks at the digital filter. It uses digital processors to carry out algorithm on different samples of the value of 
the signal. The digital filter brings arithmetical procedures on models of unmistakable time signals to 
diminish or enhance highlights of the signals. The get rid of numerous problems associated with additional 
principle class of electronic Filters such as the Analogue to Digital Converter (ADC). There are two classes 
of digital signals. That is, the finite impulse response on (FIR) and the infinite impulse response (IIR) signals. 
They are identified with the length of the impulse response. The FIR signal is an alluring choice because of 
its effortlessness in design and steadiness. By plotting the signal valves to be equivalent with the middle 
valve area, the FIR filters have a linear phase. FIR channels are recognized to contain numerous appealing 
attributes like certain consistency, the chance of exact linear phase feature at all frequencies and advanced 
execution as non-recursive developments [1]-[9]. 

Various strategies are employed for the arrangement of digital filters [10], [11]. The windowing 
strategy is the most preferred. In this strategy, the perfect impulse response is increased by being multiplied 
with a window’s function. There are various assortments of types of window capacities (Butterworth, 
Chebyshev, Kaiser and so on.). The types depend on the requirements of ripples on the pass band and stop 
band. Also, they depend on the stop band reduction and the conversion width. These distinctive windows 
bound the infinite length drive response of ideal filters into a finite window to design a genuine response. In 
any case, windowing systems don't allow enough composition of the frequency response in the diverse 
frequency bands and other filter factors such as the progress width. The Remez-exchange algorithms 
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(traditional techniques) introduced by Parks and McClellan PM does not permit express choice of the 
frequency response contained in the distinct frequency bands. Techniques based on the Remez exchange 
algorithm are the most primitive ones [12]. However, this algorithm does not give room for explicit selection 
of the top limit of the absolute ripple in the passband and stopband (õp, 5s), rather one can just determine 
their proportion. Moreover, the PM gives drifting point coefficients which requires quantization if equipment 
execution is looked for. Based on these reasons, when designing linear phase FIR filters, some digital 
stochastic worldwide model algorithm such as the Differential Evolution (DE) algorithm, Particle Swarm 
Optimization (PSO), and Genetic Algorithm (GA) are used [13]-[15]. 

The significant limitation of the designing procedure is that the similar estimations of the amplitude 
mistake in the frequency bands are specific by means of the weighting capacity, and not by the divergences 
themselves. Along these lines, if there should arise an occurrence of scheming band filters with a known stop 
band uniqueness, filter length and cut-off recurrence, the arrangement must be iterated severally. Numerous 
copies have been extended for the (FIR) signal techniques and plan designs. This is an investigative area that 
aims at accomplishing more general and spearheading strategies that are competent to determine as well as 
advance new and compound designing techniques [16]. 

Different conventional techniques exist for digital FIR channel design. Out of those one is 
the designing of the filter using the digital algorithms. This technique has a stunning capacity. That is, it 
gives alternatives of utilization of various digital algorithms and indecency in configuration. It specifically 
relies upon execution of algorithms. At first this technique was used by Parks and McClellan by utilizing 
a straightforward iterative suburbanite program and is named as PM method for filter designing. This strategy 
was altered later on by supplanting the use of basic program with enhancement algorithms. At first hereditary 
algorithm was utilized for a wide range of channel design. This was trailed by utilization of 
different algorithms. As of late, crossover algorithm and adjusted algorithms has been created from 
fundamental sort of algorithm for changes in the signal designs. Employments of enhanced molecule swarm 
advancement, versatile development molecule swarm improvements for the filter designing validate the later 
trends [16]-[24]. 

Several approaches were presented to design linear phase FIR filters. For instance, simulated 
annealing and GA method were applied to design FIR filters with coefficients values expressed as a power of 
two. However, these approaches are computationally very expensive. The Firefly method is easy to 
implement, and its convergence is controlled via few parameters. However, other methods, such as Genetic 
and DE algorithms, involves search procedures using the population genetics and natural selection process. 
The purpose of this paper is to use the Firefly approach as an alternative method to these approaches due to 
its ease of implementation in both parameter selection and the context of coding. This paper proposes a linear 
phase Low Pass FIR filter based onan alternative approach called Firefly method that assures the robustness 
of the FIR design in terms of performance and computational complexity. 

This paper is organized asfollows. In Section II, the system design and analysis of theproposed 
linear phase FIR low pass filter are presented. The Differential Evolution (DE) algorithm and the Particle 
Swarm algorithm (PSO) are discussed brieflyin Section III and Section IV, respectively. In Section V, The 
Firefly algorithm is presented. Simulation and Examples of Design FIR using Firefly algorithm presented in 
Section VI. Finally, conclusions are outlined in Section VII. 


2. THE WINDOWS FRAMEWORK OUTLINE AND EXAMINATION OF THE PROPOSED 
LINEAR PHASE FIR 
Windowing strategy is used to get the simple type of FIR filters. FIR configuration begins with a 
perfect wanted frequency response obtained as follows: 


Ha(e™) = Yin=-0 halnje“™ (1) 


Where; ha[m] is the impulse response of the distinct channel. Then, we need to make hg[m] a causal FIR 
Filterby two stages. The first hg[m] is increased by a finite length agreement d[m] (window) keeping in 
mind the end goal of getting a finite length impulse response. The second one which is causality is presented 
by time deferring the windowed drive response. One can express this in mathematical terms as follows: 


hy[m] = halm]d[m] (2) 


Where h,,[m]denotes the windowed impulse response; so, the frequency response of the windowed impulse 
response H,,(e'”) is the periodic convolution of the desired frequency response Hg(e™) with the Fourier 
transform! (e™) of the window and is given by 
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H, (e™) — = f H, (eE CEP?) do (3) 


The frequency response of the approximation filter will be 


H (e™) = e a H, (e™) (4) 


Where m; represents the necessary time delay to introduce causality of the approximated filter. Filters 
designed by windowing will have greatest error on either side of the discontinuity of the ideal frequency 


response and smaller error for frequencies away from the discontinuity. The frequency response H (e™) of 
an M length finite word-length FIR linear phase digital filter is given as: 


M—1 
H, (e™) = > hik]e™* (5) 


k=O 


Where A[K]...7k =0,1,2,...., M —1 are b-bit (sign bit included) filter coefficients. Preferable channels 


over the windowing strategy result from the minimization of maximum error yields the most preferable filters 
over the windowing method. These kinds of filters can be gotten by utilizing algorithmic procedures. In an 
effort to outline FIR filters in which some of the parameters, designed algorithms are developed. Some of 
these parameters include; the channel length (M), passband and stopband standardized frequencies (Wy, Ws). 
Researchers have developed algorithms in which M, 6,, and 6, are fixed. The other parameters have been 
enhanced. Different estimates were initially set up by Parks and McClellan in which M, wp, w,, and the 


ô 
proportion P/ Ss is fixed. From that point forward, the Parks— McClellan (PM) algorithm is the most famous 


approach for ideal FIR filters design. This is so because of its adaptability and computational productivity. 
An approximate error function in the PM algorithm is characterized by: 


EW) = GO) | H, (e™) — He”) | (6) 


Where H (e™)is the frequency response of the desired filter and Æ(e*™)isthe approximate filter, 


respectively, G (w) is a weighing function. It is used to give the weighting of the approximation error 
distinctively in various recurrence groups of frequency bands. The objective of the design is to locate the 
approximate filter coefficients that are outcomes of the ideal filter. Since the maximum error is minimized, 
then the filter is optimum. In mathematical terms, the best approximation is to be found in the following 
sense: 


(max | Ew) |) 0 


min 
{h[n]:0<n<M } 


Where A(n) is the impulse response of the approximate filter and F is the closed subset of 0 < w < 0.5. 


This model does not permit express determination of the most extreme of the outright ripples in the passband 
and stopband (6,,6,), rather one can just indicate their proportion. Moreover, the PM gives drifting point 
coefficients which require quantization if hardware implementation is sought. This motivates utilizing a 
recently developed stochastic global optimization algorithm called Firefly algorithm to design linear phase 
FIR filter in this paper. The Firefly is an insect that mostly produces short and rhythmic flashes that produced 
by a process of bioluminescence. The function of the flashing light is to attract partners (communication) or 
attract potential prey and as a protective warning toward the predator. Thus, this intensity of light is the factor 
of the other fireflies to move toward the other firefly. 


3. THE DIFFERENTIAL EVOLUTION (DE) ALGORITHM 

DE algorithm is one of the powerful evolutionary algorithms to solve real parameter optimization 
problems. Several state-of-the-art evolutionaryalgorithms use DE as evolutionary mechanism [15]-[19]. 
The DE algorithm is a population based algorithm like genetic algorithms using the similar operators; 
crossover, mutation and selection. The major contrast in building better arrangements is that genetic 
estimates depend on combination while DE depends on conversion operation. This major operation depends 
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on the distinctions of arbitrarily tested sets of methods in the populace. As a rule, DE algorithm has three 
stages [15]: 


a. Mutation: At generation k for each parent vector u, is generated according to following equation 


i rl r2 r3 
u, =x, + F(x, =x, ) (8) 


Where 5, #% #1, #1 and F is referred as the scaling factor. 
b. Crossover: This operator is essential due to its affectability to confrontdistinct and non-distinguishable 


problems. A trial child vectors yy = Ly; be y, poe ¥ : į Jare developed by the hybrid of its parent vector 


xX, and its mutant vector u, as follows: 


; E (9) 
xy if U,>C, and j#l, 


In the above model, l is anarbitrarilypicked integer from /{1; 2,...; d},U1;U2,;... ;Udare arbitrary 
independent variables uniformly disseminated in [0; 1), and C, € [0; 1] is an input parameter affecting the 
number of components to be exchanged by the crossover. 

c. Selection (Determination): A greedy Selection employed where a child substitutes its parent in next 
generation if it has a superior or equivalent fitness function value. 

In [19], they proposed the hybrid DE algorithm with local search techniques and an adaptive 
parameter value for C,. This adaptive mechanism combines the binary crossover and line recombination. And, 
refreshment mechanism is used to avoid stagnation. 

The major steps of the DE are [19]: 
1. S1: Initialize P(O) with N, individual randomly selected from the searching space. 
2. §2: Evaluate initial population P(0). 
3. repeat 
4. S3: For k=1 :MaxGen or Convergence Criterion reached 
a. Select randomly a subpopulation of N,(k) individuals with proposed solutions on P(k). 
b. Apply the fundamental DE operator “ Mutation and Cross over” to get the offspring O(k) 
c. For every child, perform a local search. 
d. Evaluate offspring O(k): if child <parent, the parent is substituted by its child 


e. If interquartile range (IQR) <V ; update the population. 
5. until stopping condition; 
6. Algorithm End 


4. THE PARTICLE SWARM OPTIMIZATION 

Particle Swarm Optimization (PSO) is an evolutionary algorithm presented and designed by 
Kennedy and Eberhart in 1995 [20]. Some efforts have been presented toward the advancement of the FIR 
Filter based on the PSO algorithm. The PSO easily applicable. Its merging might be overseen as using few 
variables. PSO is a versatile, incredible populace based stochastic pursuit or enhancement strategy with 
comprehended parallelism, which can, with no inconvenience, handle non-differential reason capacities. PSO 
is less vulnerable to get trapped on restricted optima dissimilar Genetic Algorithm (GA), replicated 
Annealing, etc. in [22-23], they expanded a similar PSO idea to present a swarm of birds. PSO is expanded 
throughout replication of bird flocking in multidimensional space. According to [24], flocking birds can be 
used to optimize a certain objective function. In PSO, every possible solution is denoted as a particle. Every 
particle is associated with a position x and a velocity v. one can express the position and velocity of the ith 
particle as following: 


atl 


=O eto) (10) 


SH 


=e Vize Vin) (11) 
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The length of every vector N denotes the dimension of the problem or number of unknown 
variables. In each iteration, the cost function is processed for every one of particles in the swarm. This 
capacity ought to be precisely designed to give a reflection of the expected outcome. The position and 
velocity of each particle are updated according to individual best, xbesti, and global best, xbest, fitness as 
shown below: 


vi” = "VE + OBL (AP — x) t epr a — x;), 


X See =x; + Atv? l (12) 


The superscripts n+1 and n indicate the time record of the current and the past cycles and B1 and B2 
are arbitrary numbers that are consistently circulated in the interim (0, 1). These random numbers are 
refreshed each time they happen. The comparative weights of the personal best position corresponding to the 
global best position are determined by the parameters c1 and c2, respectively. Both cl and c2 are typically set 
to a value of 2.0. The parameter y n is the “inertia weight” in the nth iteration. It is a number in the range (0, 
1) that selects the weight by which the particle’s current velocity relies on its previous velocity, and the 
distance between the particle’s position and its personal best and global best positions. 

The population of particles is then moved according to equation (12) and tends to cluster together 
from different directions. Since, a maximum velocity, Vmax, should not increase with any particle to stay the 
search within a desired solution space. Assessment of the cost function is done utilizing the particle’s new 
position. The algorithm goes through these procedures iteratively until the point that a particular end 
paradigm is met. 


5. THE THE FIREFLY ALGORITHM 

In [26], there is a clarification of how the Firefly algorithm that follows the firefly behavior. Firefly 
is an insect that for the most part delivers short and cadenced flashes that created by a procedure of 
bioluminescence. The capacity of the glimmering light is to pull in accomplices (correspondence) or draw in 
potential prey and as a defensive cautioning toward the predator. Along these lines, this power of light is the 
factor of other fireflies in advancing toward the other firefly. 

The light power is changed at the separation from the eyes of the onlooker. It is protected to state 
that the light power is diminished as the distance increment. The light power likewise the impact of the air 
retains by the environment, in this way the force turns out to be less engaging as the separation increment. 
Firefly algorithm originallypresented based on three idealize rules, 1) Fireflies are attracted toward each 
other’s regardless of gender. 2) The engaging quality of the fireflies is correlative with the splendor of the 
fireflies. Consequently, the less appealing firefly will push ahead to the more alluring firefly. 3) The shine of 
fireflies is relyingupon the cost function [27-30]. 


5.1. Structure of Firefly Algorithm 

In firefly algorithm, there are two essential factors, which is the light intensity force and appeal. 
Firefly is pulled in toward the other firefly that has brighter blaze than itself. The engaging quality is 
depended with the light power. The light intensity accordingly attractiveness is inversely relative with the 
particular distancer from the light source. In this manner, the light and engaging qualityis diminishing as the 
distance increment. One can express it as follows: 


2 


I(r) = I,e7Y" (13) 


I = light intensity, 
I, = light intensity at initial or original light 


intensity, 
y = the light absorption coefficient 
r = distance between firefly i and j 


Attractiveness is proportionally to the light intensity seen by another fireflies, thus attractiveness (£) is 
given by 


p=” i 
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where f represents the attractiveness at zero distance (r = 0). The distancer between two fireflies can be 
defined based on the Cartesian distance as follows: 


Lip [xi = x;| = [Bias = Xin)? (15) 

Firefly i is attracted toward the more attractive firefly j, the updated movement is expressed as follows: 
2 

Ax; = Bye (xf = xf) tae, — xf** + Ax; (16) 


where fp is for attraction, y is the limitation when the valuetends to zero or too large. If y approaching zero 
(y > 0), the attractiveness and brightness become constant, 6 = fp. In another word, a firefly can be seen in 
any position, easy to complete global search. If the y is nearing infinity or too large (y > ©), the 
attractiveness and brightness become decrease. The firefly movements become random. The implementation 
of firefly algorithm can be done in these two asymptotic behaviors. While the second the term is for 
randomization, as œ is the randomizeparameter. The £; can be replaced by ran -1/2 which is ran is random 
number generated from 0 to 1. Please refer to Algorithm 1. 


Algorithm 1:Firefly algorithm 


Input: Cost Function f (x), Initial population of Fireflies Xo, 
Define the light absorption coefficient y, Max number of 
Iteration ITR, initialize the Light Intensity Ip at xo by f (xo) 
While Loop: tillt < ITR 

For loop: for each i = 1...n all n fireflies 

Inner loop: for each j = 1...n all n fireflies 

If (I; < Ij), move firefly i towards j; end if. 
Varying attractiveness with distance r due to e~Y" 
Evaluate new solutions and update the light intensity I 

End For j 

End For i 
Rank the fireflies and find the current global best values 
End While 
Output: the best fireflies solutionsx 


6. SIMULATION RESULTS 

In this section, we show the utilization of the firefly algorithm to design optimum FIR filters for 
different cases. In all cases, for the sake of simplicity, the filter to be designed is assumed to be a linear phase 
LPF with even length to construct just 50% of the coefficients rather than all coefficients. Phase linearity of 
the filter is guaranteed by assuming symmetry of the approximate filter which reduces the dimension of the 
optimization problem to M/2. In the first case, the firefly algorithm is used to design optimum FIR filters in 
which M, wp, ws, and the ratio dp/ds are fixed. For this case, the fitness function to be minimized is 
defined as: 


E: 


cost 


max (| E(w) |). (17) 


we, UF, 


Where F, and F, are the closed subsets 0 < w < wp and w, < w < 0.5, respectively. Figure 1 presents 
the frequency response of an approximate filter with M=30, wp=0.25, ws=0.3, dp/ds=1, based on theFirefly 
algorithm, DE, PSO and PM algorithms, on the same graph. The figure shows that the Firefly response is 
precisely giving the standard FIR optimal design. This signifies that the firefly algorithm is performs well 
along the PSO and DE algorithm in terms of accuracy. In this case, the firefly, DE and PSO were randomly 
initialized and converged within a reasonable time. 
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Magnitude response of the equiripple, linear phase FIR filter 
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Figure 1. Frequency response of linear phase FIR filter for M=30, wp=0.2, ws=0.25, dp/ds=1 


For settled parameters case, the firefly algorithm is utilized to design a linear phase FIR filters in 
which M, wp, ws, dp, and ds are chosen by the designer. In Figure 2, the design process begins with the 
desired filter parameters (M, wp, ws), and the PM algorithm is employed to obtain the filter coefficients. 
Then, the designer selects a feasible value for ôp and ds based on the maximum ripple size (dpm) obtained by 


the PM.To have the capacity to control the ripples in the two bands independently, one can re-write the cost 
function in (17) as follows: 


Eos = max (| E(w) | —6,) 
È weF, 


+ max(| E0) 1-8.) (18) 


Figure 2 shows the frequency response of a linear phase FIR designed based on the firefly algorithm 
in which the filter parameters are set to M=30, wp=0.25, ws=0.3, ôp=0.1,ôs=0.01. The PM response for 
M=30, wp=0.25, ws=0.3, and ôp/ôs=100 is also shown in the same figure. Obviously, Figure 2 shows that the 
PM algorithm does not allow to control the passband and stopband ripples. However, the Firefly algorithm 
can control these qualities of ripple in the passband and stopband simultaneously. The firefly algorithm 
achieves to the desired ripples in both bands. 


Magnitude response of the equiripple, linear phase FIR filter 

E F 3 
firFirefly Design 
e  firpm Design 


m a 
© À j 
= fo ap / \s 
: He E 
= || il | 
bd Ni 
lf 
-70 
-80 
d 
-90 


0 0.05 0.1 015 02 025 03 035 04 045 0.5 
Normalized frequency 


Figure 2. Frequency response of linear phase FIR filter, forM=30, wp=0.25, ws=0.3, dp=0.1, 5s=0.01 
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It is unrealistic to acknowledge FIR channels without coefficients quantization utilizing finite word 
length. This outline requirement has not been considered in the prior work or even in the PM algorithm. 
Considering the limitations in filter implementation, the suitable word length reduce algorithm complexity 
and equipment prerequisites can prompts a dramatic change in the frequency response and design 
specification of the desired FIR. However, the capacity of firefly algorithm to deal with this sort of issues is 
exhibited in the accompanying cases. First, the impact of the quantizing process for the filter coefficients on 
the frequency response of filters designed based on the PM is demonstrated in Figure 3. In this case, the 
coefficients acquired based on the PM method are rounded by eight bits word length. As shown in Figure 3, 
the response of the PM filter with 8-bits quantized coefficients changes and its’ performance diminishes 
comparing to the floating-point PM case. In Figure 4, we present the result of the presented algorithm to 
construct the optimum FIR filter for which their coefficients are rounded using 8 bits word length for the 
same filter in previous case. 
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Figure 3. Effect of rounding on FIR filters designed by PM 


Magnitude response of the equiripple, linear phase FIR filter 


firFirefly Design 
e — rounding 


: NNN 
S sot PA S\\ 3 = ¥ 
5 |e \f \$ 1 
2 oa ; 

| | 


t t r r t r t t t 
0 0.05 0.1 0.15 02 0.25 03 0.35 04 0.45 0.5 
Normalized frequency 


Figure 4. FIR filter designed using DE with coefficients quantized using 8 bits word length for 
M=30, wp=0.25, ws=0.3, dp/ds=1 


Clearly, the Firefly performance is better than the one acquired from quantizing of the coefficients 
of the PM method. It is worth to study the impact of the order of the FIR filter on the presented method. In 
another example, the Firefly algorithm is employed to design an FIR filter with specifications M=15, 
wp=0.25, ws=0.3, dp/Ss=1 in which the coefficients are rounded using 8 bits word length. Figure 5 presents 
the performance of using the presented method to construct the optimum filter for which the coefficients are 
quantized using eight bits word length. Obviously, the Firefly method preforms well and better than the one 
acquired from rounding of the coefficients of the PM method with smaller order filter. 
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In another example, we increase the quantization bits to equal 10 bits. The firefly algorithm is 
employed to design an FIR filter with specifications M=30, wp=0.25, ws=0.3, dp/ds=1 in which the 
coefficients are rounded by 10 bits word length. The results depicted in Figure 6 show that the Firefly method 
get an enhancement on the performance and still outperforms the rounding response of PM method. 
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Figure 5. FIR filter designed using DE with Figure 6. FIR filter designed using DE with 
coefficients quantized using 8 bits word length for coefficients quantized using 10 bits word length for 
M=15, wp=0.25, ws=0.3, dp/ds=1. M=30, wp=0.25, ws=0.3, 6p/ds=1 


7. CONCLUSION 

Designing of linear phase FIR filter has been carried out through application of firefly algorithm. 
The FIR was compared with the PSO and DE algorithms. The firefly algorithm shows the ability to design 
such filters both in finite and infinite word length coefficients. There are three advantages tied to the firefly 
algorithm. That is getting to the true global minimum regardless the initial parameter, a few controlled 
parameters are used yet there is fast convergence. The results validate the fact that it is possible to design 
optimal finite word length FIR filters by use of firefly algorithm. Also, in allcases, the firefly algorithm 
performs well along with the DE and PSO algorithms. Although software implementation is important to 
investigate the capabilities of Firefly method and to simulate significant aspects of FIR Filter applications, 
hardware implementation provides real time solutions and an optimal parallelism method in terms of fast 
convergence. Hence, hardware implementations can be considered as a promising approach to implement the 
presented method which can be executed by Integrated Circuit (ICs). 
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